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Abstract 

We search for a non-SM-like CP-odd Higgs boson (a?) with m a o < in radiative decays of 
the T(IS'), using 21. 5M Y(IS') mesons directly produced in e + e~ annihilation. We investigate 
a\ — > t + t~ and a° — ► /U+jU - decay channels. No significant signal is found. We obtain upper 
limits on the product of B(T(1S) — > 7a?) and B{a\ — ► t + t~) or S(a5 — ► Our t + t~ results 

are almost two orders of magnitude more stringent than previous upper limits. Our data provide 
no evidence for a Higgs state with a mass of 214 MeV decaying to p + . Existence of such a 
state was previously proposed as an explanation for 3 S + — ► pfi + fi~ events, having / u + ' pT masses 
just above the kinematic threshold, observed by the HyperCP experiment. Our results constrain 
NMSSM models. 
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Direct searches at LEP for the Standard Model Higgs boson, a CP-even scalar, set a 
lower bound on its mass in excess of 10 2 GeV Many extensions of the Standard Model 
predict the existence of a CP-odd pseudoscalar Higgs boson (hereafter denoted as a?), which 
could be light. For example, the Next-to-Minimal Super-Symmetric Model (NMSSM) with 
a? mass below the threshold for a\ — > bb decay is particularly well motivated Radiative 
production in T(1S) decays, T(IS') — > ja®, offers a unique experimental opportunity to 
search for such a state. The couplings of the Higgs to fermions are proportional to their 
masses, therefore enhanced with respect to lighter mesons. The expected rate is given by 



where Gp is the Fermi constant, a is the fine structure constant, is the a® coupling to 
the down-type fermions, and C incorporates QCD and relativistic corrections. The coupling 
gd oc tan j3 cos 9 a, where tan (3 is the ratio of vacuum expectations for the two Higgs doublets, 
and 9a is the mixing angle between doublet and singlet CP-odd Higgs bosons (for 9a = 90°, 
a° is a pure singlet and decouples from fermions); ga depends on the detailed choice of SUSY 
parameters. For m a o < 2mb, the decay a® — > t + t~ is expected to dominate, especially at 

1 i — I 

large tan/?, B(ai — > r + r) ~ 0.9 [2j|. For m a o < 2m T , a\ — > decays are copious below 

the ss threshold. In fact, it has been suggested that 3 S + — > events observed by the 

HyperCP experiment j3] are due to a CP-odd Higgs with a mass of 214.3 ± 0.5 MeV ft 

The data for this search were acquired with the CLEO-III detector [6j operating at the 
Cornell Electron Storage Ring (CESR) and correspond to an integrated luminosity of 1.1 
fb _1 at the T(1S), yielding (21.5 ± 0.4) x 10 6 resonant decays. We also use 7 fb _1 of data 
collected at and near the T(4S) resonance for continuum background studies. 

The CLEO-III detector has a solid angle coverage of 93% of 4tt for charged and neutral 
particles. The Csl calorimeter attains 7 energy resolutions of about 2% for E 1 > 1 GeV 
and 5% at 100 MeV. The charged particle tracking system operated in a 1.5 T magnetic 
field along the beam axis and achieved a momentum (p) resolution of 0.35% at 1 GeV and 
a dE/dx resolution of 6%. The muon detector consists of the wire chambers located at 3, 
5, and 7 hadronic interaction lengths (at normal incidence) of iron absorber. 

We select events with exactly two tracks of opposite charge and at least one 7. The highest 
energy 7, which when paired with any other 7 in the event, is not within three standard 
deviations (3<r) of the 7r° mass, is selected to be a candidate for T(IS') — > •yd®. Photons 
are defined as showers that do not match charged track projections into the calorimeter. 
For optimal energy resolution and smallest backgrounds the radiative-decay 7 is required to 
be in the barrel part of the calorimeter (| cos#| < 0.8). Its lateral shower profile must be 
consistent with an isolated electromagnetic shower. When applying the 7r° veto, we place 
loose requirements on the other (non-signal) photon. Namely, the other 7 is not required to 
be isolated and the entire solid angle of the calorimeter is used for its detection. Its energy 
is required to be at least 30 MeV (60 MeV) in the barrel (endcap) part of the calorimeter. 
The 7T° veto suppresses e + e~ — > t + t~, with at least one r decaying to a channel with tt° 
(mostly pu). 

To select a\ — > t + t~ candidates we require a missing energy between 2 and 7 GeV. The 
total energy calculation is based on charged track momenta (the pion mass is assumed) and 
calorimeter energy for the looser definition of the 7 candidates (see above). To suppress 
hadronic events from the continuum production and Y(IS') decays, at least one of the 
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charged tracks must be identified as an electron or a muon. Events with two electrons are 
discarded to suppress Bhabha like events. Besides being identified in the muon system, the 
(jl candidate is required to have an energy deposited in the calorimeter (E) consistent with 
a minimum ionizing particle. The electron candidate must have E equal p within ±15%. 
dE/dx consistency is required for for leptons. The invariant mass of photons (except for 
the radiative-decay 7) plus the charged track not identified as a lepton is required to be less 
than 2 GeV. To suppress final state radiation, the cosine of the angle between any charged 
track and the 7 candidate must be less than 0.99. 

We search for evidence for a signal of a monochromatic peak in the 7 energy distribu- 
tion. Thus, our results assume that the a\ natural width is negligible compared with the 
experimental resolution, an assumption which is expected to be true, with the exception 
of the heaviest masses probed in the t + t~ channel, due to possible mixing with the rjb Q. 
CLEO previously published an alternative method for probing the a? mass approaching the 
bb states, which is not sensitive to assumptions about its width [§]. 

The selected event sample is composed mostly of continuum e + e~ — > (7)t + t~ events, 
where the 7 candidate comes either from initial state radiation (ISR) or from a 7r° produced 
in r decay, with the second 7 not reconstructed. The background estimates are superim- 
posed on top of the spectrum obtained at the T(IS') resonance in Fig. QJ1. The continuum 
backgrounds are estimated by scaling the T(4f>) distributions. There is also a significant 
contribution from T(1S) — > t + t~ with the 7 candidate coming from a ir° decay. The ob- 
served 7 spectrum with binning comparable to our 7 energy resolution is shown in Fig. [lb. 
No significant peaks are observed. 

The channel a\ — > is selected by identifying both muons. We require that the total 

observed energy be within 250 MeV of the center-of-mass energy. The invariant dimuon 
mass has better resolution than the 7 energy, therefore we use it to look for the a\ signal. 
The selected data are dominated by radiative /i-pairs with a hard radiative photon. The 
data selected at the peak of the T(IS') resonance are well described by scaling the data 
collected at and near the T{AS) as illustrated in Fig. [lb. The T(15) distribution plotted 
with binning comparable to our dimuon mass resolution is shown in Fig. [TH. No significant 
peak is found except for the J/ip produced by ISR. 

The signal efficiency varies with Higgs mass, or equivalently, 7 energy. In order to de- 
termine the efficiency, we generated signal Monte Carlo (MC) for several a\ masses and 
interpolated for masses in between. Proper angular correlations were implemented in the 
MC for both the polar angle of the radiative 7 and for r polarizations [9j . The signal peaks 
observed in the dimuon mass (7 energy) distribution for a\ — > (a§ — > r + r~) were 

fitted to a Gaussian (with an asymmetric low energy tail, i.e., a Crystal Ball line shape [ToT ]) 
to determine reconstruction efficiency and detector resolution, which are shown in Fig. [Sb.c 
(bd). The fit to the J/ip — > ISR peak observed in the T(IS') data gives a value for 

the resolution consistent with the MC expectations for an a\ — > ^ \i~ signal at that mass 
(Fig. HI). 

We have scanned the observed 7 energy and dimuon mass distributions by fitting a signal 
peak on top of a cubic background polynomial, changing the peak position in steps equal to 
the detector resolution. The peak width was fixed to the MC expectations. The fit range 
was set to ±0.5 of ln(E 7 ) (±0.25 of ln(M w )) around the peak position. The J/ip — > 
peak region was excluded within ±3cx, unless fitting a peak at the J ftp mass. Since, in 
the dimuon channel, the continuum backgrounds saturate the T(1S) sample (see Fig. [lb), 
we simultaneously fit the background polynomial to the dimuon mass distribution obtained 
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from the higher statistics T(4S) data. The overall background normalization factor for the 
T(4S) data was included as a free parameter in the fit. 

To test for possible bias in the fit procedure, we calculated average and root mean square 
(RMS) values of fitted signal amplitude (N) divided by its error (AN). In absence of any 
signal peaks, values of 0.0 and 1.0 are expected, respectively. The average value for the data 
is +0.01 ± 0.09 (-0.06 ± 0.05) for the 7t + t~ (7//+^") sample, while the RMS is 1.16 ± 0.09 
(1.05 ± 0.05). To cover the observed deviations from the expectations we assume 15% for a 
possible systematic error in the fit procedure. 

The ratio of the fit likelihoods for a signal peak included in the fit (C max ) and the data fit 
with the background term alone (£0) is used to calculate the peak significance in standard 

deviations, \jl \n(£ max / £ ). No peak with significance above 3a is found in the •yr + r~ data. 
In the 7/x + /i~ data, the ISR J/if> peak has 8.3a significance. There are two other mass points 
which produce peaks with significance slightly above 3a: 3.3a at 2041 ±4 MeV and 3.1er at 
211.92 ± 0.15 MeV. For one trial, the probability (e) of background fluctuations producing 
a peak with significance equal or larger to 3.3a (3.1a) is 0.05% (0.1%). We performed 482 
fits to 7/i + /i~ spectrum with the peak positions separated by one unit of mass resolution. 
Assuming that peaks must be separated by at least 3 units of mass resolution to fluctuate 
independently, we performed about n = 482/3 ~ 160 statistically independent trials. The 
overall probability in our scan of producing at least one peak with significance of at least 
3.3a (3.1a) is 1 - (1 - e) n « 8% (15%). 

With no evidence for an a® signal we set upper limits on its possible production rate 
as a function of the a® mass. To determine upper limits on the signal event yield we 
fix the signal amplitude to positive values and minimize it with respect to the background 
parameters to obtain a likelihood for given signal amplitude. We then integrate the likelihood 
function and find a signal amplitude which bounds 90% of the total area. Dividing this event 
limit by efficiency and the number of T(1S) decays in our data we obtain upper limits on 
B(T(1S) — > 70?) x B(a\ — > l + l~) (I = r or //), which are displayed in Fig. [31 The limits 
were scaled up by 20% to account for the possible systematic error in the fit procedure 
(15%), in the efficiency calculation (MC statistics, interpolation between MC points, detector 
modeling; together < 10%), in number of T(IS') decays (2%), and in simulation of detector 
resolution (10%). Our upper limits on B(T(1S) — > ja®) x B(a\ — > r + r~) are almost 2 orders 
of magnitude more stringent than previously obtained by ARGUS [ll[ . Our upper limits on 
B(T(1S) — ► 70.?) x B(a\ — > yU + /i~) are the first experimental bounds. 

Our 7t + t~ results provide stringent constraints on NMSSM models, eliminating a large 
portion of previously unconstrained parameter space. This is illustrated in Fig. HI where 
NMSSM calculations by Dermisek, Gunion, McElrath [2} are compared to our upper limits. 
While some model parameters were fixed (e.g. tan/3 = 10) in the theoretical calculations, 
other parameters were sampled. Each point represents a different choice of NMSSM param- 
eter values consistent with the current experimental constraints. The plot on the right rep- 
resents models with the additional requirement of low fine-tuning of electroweak symmetry 
breaking (see Ref.jl] for details). Color coding corresponds to various a\ mass ranges. Our 
upper limits in various a\ mass ranges are shown by horizontal lines. Solid (dashed) line(s) 
represent an average (minimal and maximal) upper limits in given mass range. Assuming 
that B(a° 1 — > t + t ) is 100%, the models above these lines are excluded. Our upper limits 
in the 7t + t~ channel for lower a® masses, 2m T < m a o < 7.5 GeV (red lines at the center of 
each plot), eliminate a significant fraction of models in this mass range (red points). Only 
very few models are challenged in the 7.5 < m a o < 8.8 GeV range (green points and lines at 
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the sides of scatter plots). For higher masses (black points and lines touching the plot axes) 
our discriminating power fades away as the backgrounds increase while the expected signal 
rate decreases due to the phase space suppression. The non-singlet fraction of a? (cos #,4) 
increases with falling tan/3 0], though the net effect on g d , and therefore B(T(1S) — > 70°), 
is to decrease. At the same time, a° coupling to up-type fermions, g u oc cos 9 a/ tan (3, in- 
creases, lowering B(a\ — > r + r~). Thus, models with small tan/5 values are less constrained 
by our data. 

Since B{a\ — > yU + /i~) is not expected to be large above the a\ — > ss threshold (~ 1 
GeV), we did not transfer our limits on B(T(1S) — > 7a°)B(a° — * /i + to Fig. HI where 
NSSM model calculations of B(T(1S) — > 70°) below the r + r _ threshold (blue points) were 
performed only for m a o > 1 GeV. Our limits below this mass value constrain NSSM scenar- 
ions. Of particular interest is our upper limit for m a o = 214.3 MeV, i.e., the /i + yU~ mass of 3 
E + — ► j9/i + yU~ events observed by the HyperCP experiment The fit to our data (see Fig. [5]) 
gives 7.0^^5 events at this mass and leads to an upper limit of B(T(1S) — ► 70°) < 2.3 x 10~ 6 
at 90% C.L. He, Tandean and Valencia showed that they could explain the HyperCP events 
with the a? hypothesis and still be consistent with the constraints from K — > 7Tfi + fi~ ex- 
periments [5|. In their calculations they used g\ = 0.12, while our upper limit translates 
into g\ < 0.026 (using C = 0.5 [3]), which calls for reevaluation of the a? hypothesis for the 
HyperCP events. 

In summary, we have obtained meaningful upper limits on B(T(1S) — > 70^) x B{a\ — > 
t + t~) and B(T(1S) — > 70°) x B(a!{ — > ji^). Our limits on 7r + r _ are almost two orders 
of magnitude more stringent than those from ARGUS and eliminate a large portion of pre- 
viously unconstrained parameter space in NMSSM models. Our limits on 7/i + yU - challenge 
models with a\ mass below ss threshold and the a\ interpretation of HyperCP S + — > pfi + fi~ 
events. Our limits are applicable to any light scalar or pseduo-scalar boson, which arises in 
various extensions of Standard Model. 
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ties Council. We thank Radovan Dermisek, Jack Gunion, Michelangelo Mangano, German 
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FIG. 1: Photon energy and dimuon mass distributions in 7t + t~ (a,c) and 7/i + ^ _ (b,d) data, 
respectively. Bin size for the right column plots is given in the axes labels in parentheses. In 
the top row, the T(IS') data (points with error bars) are compared to the estimated backgrounds 
(dashed and solid lines). In the bottom row, the T(1S) data (solid line) are shown in fine binning 
comparable to the detector resolution (see bottom row of Fig. [2]). In (b) the J/tp ISR peak is 
shifted in the background estimate, since we scaled /x momenta down by the ratio of the beam 
energies when scaling the higher energy data to the T(15) distribution. 
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FIG. 2: Efficiency (a,b) and a\ mass resolution (c,d) obtained from the fits to the a? — * t + t~ 
(left column) and a? — > fJ^ + fJ-~ (right column) signal MC (points) and interpolated for the regions 
in between (solid line). In (d) relative dimuon mass resolution was multiplied by a factor of 10. 
See Appendix C of Ref. Q for explanation of improvement of the dimuon mass resolution near 
the kinematic threshold. The hollow point with the error bar in (d) represents the fit of the mass 
resolution to the J/tp — > /i + ISR peak observed in the T(IS') data. 
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FIG. 3: Upper limits on B(T(1S) 7a?) (a) x8(a? -» fj, + /i~) (b) xB(o? -» r + r~) as a function 
of the a® mass (90% C.L.). The color coding corresponds to the one used in Fig. SJ The dashed line 
indicates the region (m a o > 9.2 GeV) where a® is likely to mix with 775 and acquire a non-negligible 
width, thus invalidating our analysis method. 
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FIG. 4: Comparison of CLEO upper limits on B(T(1S) — > 70°) x B(a\ — > r + r~) (solid and dashed 
lines) to the NMSSM predictions by Dermisek, Gunion, McElrath (points) [2(. See the text for 
explanations. 
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FIG. 5: A fit of a peak at a dimuon mass of 214.3 MeV with fixed width at the expected mass 
resolution, on top of cubic polynomial to our 7/i + /i~ data obtained at the T(15) center-of-mass 
energy (a). The polynomial describing the continuum backgrounds was simultaneously constrained 
by the data collected at and near the T(45) resonance (b). 
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